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The paper deals with various approximate procedures for calculating the distortion of a pulse 
after it has propagated through a dispersive channel such as a waveguide. The methods used for 
evaluating the integrals utilize a stationary phase principle. Both impulsive-type and quasi-monochro- 
matic sources are considered. It is shown that, in most cases, the transient response may be obtained 
straightforwardly from the shape of the phase versus frequency characteristics of the system. Some 
attention is given to the complication which arises when the group velocity has an extremum as a 
function of frequency. 

1 . Statement of Problem and. General Introduction 

The problem is to calculate the waveform of an electromagnetic pulse after it has propagated 
through a linear dispersive channel or medium. The waveform of a component of the electro- 
magnetic field at the input is designated Eoit) and the transformed waveform at the output is 
E{t). Thus, E(){t) is given, while E{t) is sought. 

To facilitate the solution, the source field is written as a Fourier integral 

Eo{t) = ^ r\{o))e^-Ui(o, (1) 

^rr J -^ 

where, for the moment, g{oj) has no singularities on the real axis of the complex co plane. The 
inverse of (1) is given by 

g{oj)=r E,{t)€-'-'dt, (2) 

J —oc 

where g (co) may be regarded as the spectrum of the source pulse. 

After the spectral component g((jt)) exp {iojt) passes through the system it will be modified by 
the transfer function Rico) exp [—iip{aj)] where R{o)) contains the amplitude and any other slowly 
varying complex functions of co. Thus, the output pulse is represented by 

1 r+^ 

^W = ^7~ RMgioj)e'^^e-'^^'^^d(jo. (3) 

Ztt j -X 

In some cases, the evaluation of this integral may be carried out by analytical means when the 
functions gico), R{oj), and (^(co) have simple algebraic forms. Unfortunately, more than not, the 
transfer function R(q)) exp [— iip{o))] is some comphcated function of frequency which may be 
given only in numerical form. In such cases, it may be feasible to evaluate the integral by a 
purely numerical procedure. With the wide availability of the digital computer, this is certainly 
fashionable at the moment. Consequently, one might say that the problem has been solved and 
no further discussion is needed. However, it would be a pity if one accepted this answer since 
all physical insight into the nature of transient processes has been ignored. Furthermore, apart 
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from the economic considerations of using large-scale computers, the direct numerical procedures 
require very special precautions when the integrand of (3) is highly oscillating. The awkwardness 
of the situation is compounded when the numerical data for R{aj) and (p((x)) are given only at a 
finite number of real frequencies, let alone the possibility that data for a whole frequency range 
might be missing. 

From the above considerations, there seems to be a need to utilize approximate procedures 
for handhng integrals of the type given by (3). Furthermore, by exploiting some of the rapidly 
varying characteristics of the transfer function, certain approximations can be made which lead 
to great simplifications. 

The point of view adopted here is similar, in many aspects, to those expounded many years 
ago by Sommerfeld [1914] and Brillouin [1914]. However, in this paper, the emphasis will be on 
transfer functions which are characteristics of propagation channels such as the earth-ionosphere 
waveguide. Also, we shall be concerned with source pulses which have spectra covering a wide 
range of frequencies. 

2. The Stationary Phase Method 

As a starting point, we write the integral (3) in the following form: 

y= I G(co) exp [—iF{a))]da), (4) 

J— 00 

where Fico) is real and contains all the rapidly varying phase terms, while G{a)) is a slowly varying 
function which may be complex. In accord with the classical notions of the stationary phase 
method, the important contributions to the integrand are where the function F{(o) is stationary. 
This suggests that we write F(co) as a Taylor expansion about the frequency (Ds which is defined by 

F'{cos)=[dF{co)ld(o]^=^=0. 
Thus, 

F{o}) = F{a)s)^ ^ F(w,)H 7— ^F'"(w,)+. ... (5) 

For many purposes, the series may be truncated beyond the term containing the second derivative 
F"{a>s). Furthermore, in the physical systems under consideration, the stationary points usually 
occur in pairs at co = ± ojs. 

The contribution from + cos may be approximated in the following manner: 



J^=G{o)s)e-i'^^-^s'> j^^exp [-i{(0- (OsfF'MI2]da), 



(6) 



where the slowly varying G(co) in (4) is taken outside the integral and replaced by its value at co = cos. 
The integral is now easily transformed to the well-known form 



J -00 



e-^'dx^MiylK (7) 

J — 00 

Thus, 



J+ = G{(Os) 7^ 



iF'iojs) 
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1/2 

exp [-iF{(Os)l (8a) 



In a similar manner, the contribution from oj = — a)s is found to be 



J- = G*(co,.) 



277 



1/2 



where the asterisk denotes a complex conjugate. Writing, 

G{ojs)= \G{cjOs)\ exp [//?(co.s)|, 



exp f-h/F(co.s.)], 



it is a simple thing to show the combined contribution is 

(277)1/2 

J = J+-^J- = 2\G{cOs)\ T-r^ — -r— COS 



77 



F(0)s)-g[(x)s)±- 



(8b) 



(9) 



(10) 



where the -^rrl^ is to be used when ficos) > and the —77/4 is to be used when f"(a>,s) < 0. 

When the integral has additional saddle points, the c()ntril)ution from them must be added to 
those from ± cjs- 



3. Pulse Propagation in an Ideal Waveguide 

The applicabihty of the simple stationary phase method is best illustrated by an example. 
The model we shall choose is a vertical electric dipole sitting on the bottom surface of a parallel 
plate waveguide. In terms of cyhndrical coordinates (p, (/), z), the dipole of length ds is located 
at the origin, while the perfectly conducting walls are dehned by 2 = and z = h. When the dipole 
current varies as /(w) exp {icot), the magnetic field has only a (j> component given by [Wait, 1962] 



^^(^)^~^j^ — ' ^ S e,MUs„f>) cos ac„z) 

^^ ^P „=o, 1,2, .. . 



(11) 



where eo = 1 , en = 2(n t^ 0), 



Sn = il-Ciyl^ Cn = 7Tnl{kh), 



and k = (o/c. The function H^f^^ is the Hankel function of the second kind of order zero. Now, in 
this particular problem, the input function is the waveform ioit) of the current in the source dipole 
and the output is the magnetic field h^(t) as a function of time t. These are related to the time- 
harmonic quantities /(co) and Hd^ioj) in an analogous manner to (2) and (3). Specifically, 



ct>(co)e''^^d(o 



and 



1 r+^ 

I{(o)= I " io{t)€-'''^dt. 

J —oo 



(12) 



(13) 



For this example, we shall choose ioit) to be a suddenly applied impulse (e.g., an idealized lightning 
discharge). Thus, 
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(14) 



where d{t) is the unit impulse function. Thus, from (13), it follows that: 



(15) 



which, in effect, says that the source has a uniform spectrum over all frequencies. 

The transient response of the waveguide channel for this impulsive current may thus be written 



where 



lods d ^ /7Tnz\ 



^^•"-iI>-''4h'-(T)H 



(16) 



d(o. 



(17) 



The integral Vnip, t) will now be evaluated by the stationary phase method and the result 
will be compared with the known exact result. The central assumption is that the Hankel function 
may be replaced by the first term of its asymptotic expansion. Thus, 



Vu{p. t) 



2J-00 



C(a))e-^^(-)rfco, 



where 



and 



G((^y- 



2A"' 






-1/4 



F(oy) 



77/7 C 



1/2 , 



- (x)t. 



(18) 



(19) 



(20) 



The radical in the latter expression is to be chosen so that the real part is positive if | w | > | TTncjh \ 
and the imaginary part is negative if | co | < | Trnclh |. Physically, these conditions amount to say- 
ing that only outgoing propagating modes are permitted for p-^°<^, while the cutoff modes are 
damped in the positive p-direction. 

The stationary phase point, obtained by applying F'{(x)s) = to (20), is given by 



where 



cos = {7rnclh)ctlX 



X=[(c^)2_p2]l/2, 



It is also a simple matter to show that 

FM = - irnX/h, F\ws) = - hX^KTrncPc^), 
and, finally, 
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(21) 



where u{t')= 1 for ^' > 0,==0 for t' < 0, and t' = t — {p/c). This is to be compared with the exact 
value given by [McLachlan and Humbert, 1941; Budden, 1951] 



2c 



7TnX\ 2 . (TTnX\ , ct-\rX 



u{t'). (22) 



When X<< ct or p, and if p > > /?, the term containing the logarithm is negligible. On close 
inspection of the matter, it is seen that these are just the conditions for the stationary phase method 
to be valid. It is also worth noting that (21), for h = oo^ may be checked by comparing it with pair 
912.3 in Campbell's and Foster's [1948] tables. 

The transient magnetic field, within the stationary phase approximation, is found from (16) 
and (21) to be 






/l=0 



COS 



7rnX\ , f7TnX\ . /7rriX\ 1^^ (irnzX . 



V /, / ■■'■ ' (23) 



which is valid when X «ct or p and p» h. 

It is clear from (23) that the individual modes, for /i > 0, vary with time in a quasi-sinusoidal 
manner with a quasi half-period which is approximately given by 

T U 1 

This shows that the apparent period of the oscillation increases slowly with time. It is also in- 
teresting to note that the higher order modes have shorter oscillation periods. Of course, in the 
case of the zero-order mode (i.e., ai = 0), the magnetic field varies with time monotonically as 
[(c^)2-p2]-3/2. 



4. Application to More Realistic Models 

As indicated above, the application of the stationary phase method is quite straightforward 
when apphed to transient calculations in conventional waveguides with no wall losses. On the 
other hand, if the waveguide walls have finite conductivity, the problem may be very difficult if 
highly precise results are desired. Fortunately, for most cases of practical interest, the attenua- 
tion of the modes is accounted for by a multiplicative factor which, in the frequency plane, may be 
designated by exp [— Q:/t(co)p] where an{o)) is the attenuation coefficient as a function of frequency. 
Within the approximations of the stationary phase approach, the transient response of the indi- 
vidual modes is then simply multiphed by exp [— a;j(a>s)p] where cos is the stationary phase point 
as computed for the lossless case. This perturbation technique has been applied on earlier 
occasions to the computation of transient waveforms in the earth-ionosphere waveguide [Wait, 
1958]. Apart from the plausibility of the approach for small losses, a mathematical justification 
has been given by Gajewski [1958] for metallic waUed waveguides at centimetric wavelengths. 

In order to discuss the full significance of the stationary phase method in mode theory, it is 
desirable to return to (4). When the integrand represents an individual mode, the rapidly varying 
phase function F[(i>) may be written 

F(a>)-^-co^ (25) 

under the assumption that the phase characteristics of the source and other factors are sufficiently 
slowly varying to be lumped into G(co). Here, it is understood that Vpi^) is the phase velocity of 
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, given mode. The stationary phase condition may now be written 



F'M--^-^ = 0, (26) 



is, by definition, the group velocity. Often it is convenient U 



where Vgico,,) = 1 / ,. . . 

use (26) to estimate the stationary point cos from given data on fy(a)) plotted as a function of o). In 
a similar manner, one may wish to estimate the second derivative F"{cos) which occurs in the 
denominator of (10). In this case, 



FX(Os) = p 



1 



d(x) Vg((0) 



(27) 



which, as indicated, is the derivative evaluated at co = (Os. 



In the preceding development, it has been tacitly assumed that the phase function F(co) may 
be adequately approximated by the first three terms in a Taylor expansion of F{oj). Clearly, 
this is inadequate when the second derivative F"{(Os) vanishes. This would lead to particular 
difficulty in trying to evaluate (10). As an alternative for situations of this type, we develop F(a>) 
in an expansion about the point ojg which is defined by the solution 

F"i(o) = ^^F{co)=0 (28) 

or 

d I _ 1 dvgia)) 



d(x)Vy{o)) Vg(o)) d(o 



0. (29) 



This condition corresponds to the situation where the group velocity, as a function of frequency, 
has a minimum or a maximum. In the case of practical waveguides in radio technology, Vg{(o) is 
usually a smooth monotonic function and the condition is not satisfied for any finite frequency. 
However, there are numerous examples in seismology [Ewing, et al., 1957] and, in certain cases, in 
terrestrial radio waveguides [Wait, 1962] where the group velocity has a maximum or minimum. 
For example, an analytical model of the earth-ionosphere waveguide would indicate the group 
velocity, for the dominant mode, has a maximum at a frequency (wg/27r) of about 22 kc/s for typical 
daytime conditions [Wait and Spies, 1964] (i.e., where the phase velocity versus frequency curve 
has maximum curvature). 

Instead of (5), we now write 

F{co) = F{cog)^iw-ojg)F'{a)g)^{o)-a)grF'"{a)g)l6^. . . . , (30) 

where the second derivative term is absent on account of 

F"{<Og} = Q. 

When the series in (30) is now truncated beyond the third term, it is a simple matter to show that 

y+ = 2G(a)3)e^'«'"9» f " cos (d.F ' + ^^) dO.. (31) 
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^here 



il = oj — ojf,. F' = F'io),,), and F'" "= F"'(ci),/). 
When F"' > 0, the approximate expression for J+ may be written in the following form: 

/ 2 vi^ r / 2 V^l 

A^277i/2GK)(^— j e-^^^^-yh^F' i^—j J' 

where v{Z) is the Airy integral defined by 

v{Z) = — [ " cos [Z5 + (s^l3)] ds. 



In a similar manner, we find that the contribution from a) = — co^; is given by 



/ 2 \^l^ 
J- = 27ri/2G*(co,) f ^j e'-'^'Wy 



■W] 



As a result, 



J = J++J~ = — rpTn^ — COS [/'(co^j-^wyjji; 



1 p,, 



where, in analogy to (9), 

G{(o,j)=\G{a)ff)\ exp [i^cov)]. 
A very similar expression holds for F"' < 0. Then, it is convenient to write 

/ 2 \^l''^ 






When this is applied to a waveguide mode, 



-F' 



.f" 






(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



where (Og is defined as the solution of 



Also, for this situation, 



and 



do) 



F' 



Vg{(x)) = 0. 



Vg{0)c 



-—t. 
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Figure 1. The Airy function and a related integral. 
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From the above, it is evident that the transient response of a waveguide mode is an oscillatory 
function with a constant frequency co^/27r. It has an amplitude which is proportional to |F'"|~^/^ 
and the envelope is proportional to the Airy function v{ — x) where x = — F\2lF"'yl^. It is under- 
stood that the transient response so calculated is valid only when x has a magnitude not large 
compared with unity. If this is not the case, higher order terms in the expansion in (30) should be 
considered. To show the manner in which the envelope of the pulse varies with time in this region, 
the function v{ — x) is plotted in figure 1 as a function of x from — 5 to + 3. For a waveguide mode, 



x = 



Vgi^^g). 



2 y/^ 

F'") 



(38) 



which is positive when t > p/vg{(Og) since F'" is positive when the group velocity has a maximum 
(at a) = ojg). Thus, in a sense, the curve of v( — x) versus x in figure 1 may be regarded as the time 
history of the envelope of the leading portion of the pulse response of the waveguide. In this case, 
time increases from left to right. On the other hand, if the group velocity has a minimum, F'" is 
negative and then we should note that 



VgioJy), 



-f" 



1/3 



(39) 



which indicates that time increases from right to left. 



5. The Quasi-Monochromatic Pulse 

The approaches mentioned above are applicable to propagation of the broadband pulses in a 
dispersive channel. In these cases, the spectral characteristics of the source are sufficiently broad 
that they may be lumped into G(co) and taken outside the integral in either (4) or (31). However, 
when the source has most of its energy in a relatively narrow band of frequencies, another approach 
is desirable. 

An important practical situation corresponds to choosing the source to be of the quasi-mono- 
chromatic form 



E^{t) = A{t)e^^^\ 



(40) 



where coo is the (angular) carrier frequency and A{t) is a specified well-behaved function which may 
be called the modulation envelope. The corresponding spectrum, using (2), is 



^w) = A{t)e'^'^o-^)tdt. 
1394 



(41) 



Using (3), the response of the dispersive propagation channel is then to be obtained from 

Fit) = ^-1^1^^ A(t)e^^--^^^^^'-^^^^^R{(S)do}di, (42) 

where /?(a>)e~**^^^ is the complex transfer function of the channel. The phase function (^(co) is now 
expanded in a series about the carrier frequency ojo- Thus, 



(^(w) = (^(w,)) + IV(wo) + (n2/2)(^"(coo) + (a3/6)(^'"(a)o), 



(43) 



where fl = aj — coo. It is tacitly assumed that higher order terms in fl may be neglected.^ 

We now seek a method to carry out the w integration in (42). To facilitate this step, it is 
convenient to replace exp [— i(n^/6)(/:>"] by a suitable integral representation which involves only 
the first power in H in the exponent. This is obtained by noting that the definition of the Airy 
function v{y) of argument y may be written as a Fourier transform [Wait, 1962] 



2ttJ- 



e-"'l^e-"<idx =—rjr, v{y). 



The inverse of this transform is 



^ir^ 



,(y)f.'>.Vjy=e-'-^"'/'\ 



Using this result, it is then readily found that 



exp 



4^ r\iy) exp [i{^'''l2yf''yn]dy. 

** J — oc 



which is the desired form. 

On using (43) and (46), it is now seen that (42) takes the rather ominous-looking form 



£'(0 = e'^<''" 



Mw()) • 



1 



We now introduce a new variable z such that 



n+ ^"-^+^'-!^"'^^)'"^ )(^T^--^-/-.. 



The exponential term in the integrand of (47) then becomes 



exp 



■i('!TZ^I2) + i 



{i-t+<p'-{<p"'i2yi^Yf 

2(p" 



(44) 



(45) 



(46) 



r+oc r+oc r+oc 
2 A(t)v{y) 

J — oc J —oc J —X 

X exp iO.[t-t-<p' + {<p"'l2yl^y]-^ jflV'l dydida. (47) 



(48) 



(49) 



Now, since we are dealing with frequencies co near coo, the hmits of the new variable z may be re- 



' It is also assumed here that /?(a>) is slowly varying compared with exp [— i<p{(i})] in the vicinity of wo. Thus, it may he taken out of the intejirai and replaced 
by RioJi)). In most of what follows, we shall set R((iJo)= 1 for the sake of convenience. 
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garded as — ^ and + oo. Thus, using the resuh that 

exp [i7TzV2]dz = 21/2 exp (^77/4), (50) 



i: 



we obtain 

£(0 = e'W-Ma>o) —L— r^ f ^^ Ait) exp [^(f- ^ + (^^ - ((^^ " l2y/'yYI{2^")] X J nr) JyJ^\ (51) 

i2iyi^iT j-x J-x V(^"(w()) 

Then, on introducing another variable u defined by 

Ct-t+<p'-{<p"'l2yl^yf _^ , 
2^" 2 " ' 

it is seen that (51) becomes 
E{t) = exp [i(w()f — (p(o)(,) — (77/4))] 

~-J^^ i^''v(y)A[t-<p' + {<p"'l2yl^y+{7T<p"Yl^u] exp {inu^mdudy. (52) 



^(2. 



This general expression, although it appears to be rather complicated, may be used as the starting 
point for a number of investigations dealing with the effect of dispersion of pulses. For example, 
if it is permissible to set ip"' = 0, (52) reduces to 

1 r+^ 

E{t) = -a-i)e^^of-i^^^o^ A[t-ip'-\-i7npyl^u] exp {i7Tirl2)du (53) 

■^ J —X 

by virtue of the identity 

r+x 

I v(y)dy=-7T^I^\ (54) 

which is a special case of (45) when x = 0. Equation (53) is identical to the special case derived 
by Ginzburg [1964] who neglects cp"' at the outset. If, in addition, (p" = 0, it is seen that 

Eit) = e'^o^-Ma>o) A[t -cp'l (55) 

which is the well-known result for a nondispersive channel. Here, the pulse envelope has been 
delayed by (p'icoo) seconds but has not changed its shape. On the other hand, the carrier phase 
is determined by ipioJo). The quantities (f' and (f are referred to often as the group delay and the 
phase delay, respectively. 

6. Discussion of Pulse Distortion 

Equation (53) may be used to study the distortion of the pulse shape, provided the (f'" term 
is negligible. For example, if the original pulse has a rectangular shape such that 

A{t) = OU)T t<-Tol2 

= 1 for t > - Tol2 and t < Tol2 

= 0{ort>TJ2, (56) 

1396 



it readily follows that: 



m=^e> 



fu., 

J Hi 



exp {iTnrl2)dn, 



(57) 



where 



U\- 



-e 



[V'(wo)]^/^^ 



Uz- 



To-0 



and 6 = {To/2) -\r t — if' {(Oo). Here, 6 is the time measured from the instant — (Tol2)-\- (p'(a)o). If the 
dispersion is sufficiently small, Ui and uz become infinite and (57) reduces to (55), as it must. 

To illustrate the influence of the dispersion, it is assumed that the sitj;nal duration To is large 
compared with [7Tip"{oj{))y^^. Then, the form of the leading edge of the signal may he written 



\m\ 



1 



/■ 

J III 



exp {imrl2)(lu 



(58) 



where Wi = — 6^/|7r(/?"(cO())r^^. Using readily available tabulated data for the above Fresnel integral, 
E(t) is plotted as a function of— ii\ in hgure 2. The abscissa, in this case, may be regarded as a time 
scale where zero corresponds to the time t = ip'{o)()) which is the total group delay. As indicated, 
some energy arrives before this time for an actual dispersive medium. The time for the pulse to 
actually l)uild up to its final value may be described conveniently by a parameter ti,. This is defined 
as the time from ^=0 for the pulse envelope |A>(/)| to a{)[)roach within 5 percent of unity. From 
inspection of figure 2, it is seen that 



th = 4[7r(/?"(aj())p^'- seconds. 



(59) 



A tyi)ical value of //> for daytime propagation in the earth-ionosphere waveguide of height 70 km, 
at 15 kc/s, for a range of 5000 km is .']()() /xsec. The parameter' //; is proportional to the square root 
of range and it depends in a rather complicated maimer on other geometrical |)aranieters. How- 
ever, it is only slightly dependent on the nature of the ionospheric model assumed once the reflect- 
ing height has been specified. 

It is interesting to note that the distortion of the envelope shape indicated in figure 2 is very 
similar to the spatial pattern resulting from diffraction at a semi-inhnite screen. In general, for 
the waveguide case, decreasing the carrier frequency towards cutoff causes the envelope function 
to spread out increasingly (i.e., th becomes larger). Also, as may be seen from a detailed study 
of the complex form given by (56), the instantaneous frequency shows a modification from the 
carrier coo- This may be described as "angle modulation." 



Figure 2. The leading edge of the envelope of a step- 
modulated carrier in a dispersive medium. 




-2-10123 
Time Parameter, 9/y¥^yIj^ 
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Using essentially the present second-order theory (i.e., neglecting <^'", etc.), a number of 
papers have appeared which deal with distortion of pulses both in acoustic [Pearson, 1953; Proud 
et al., 1956] and electromagnetic waveguides [Elliott, 1957; Karbowiak, 1957; Kotvun, 1958]. A 
rather disconcerting aspect of Elliott's [1957] work is that his calculated envelope patterns do not 
have the oscillations near the final buildup region of the pulse. The existence of these "Fresnel 
wiggles" has been recently confirmed by Ito [1964] who carried out some very convincing experi- 
ments in an electromagnetic microwave guide at frequencies of 7.8 and 8.2 Gc/s, which were 
respectively 1.19 fc and 1.25 /f where /c is the cutoff frequency. Similar effects were observed 
by Walther [1961] who worked with acoustic waves in a water-filled channel. 

7. Complications at an Extremum of the Group Velocity 

The generalization of the above treatment for the rectangular pulse to allow for the third deriva- 
tive (p"'{o)o) may be readily obtained from (52). Thus, using the special form of the source envelope 
given by (56), it is a straightforward matter to show that 



E{t)-- 



1 



27r>/2 



_^ ^i(DQt-i<fi(t)Q) 



1: 



^here 



To-d-{ip"'l2yl^y 



(V)l/2 



-fF 






v{y)dy. 



(60) 



F{ 



jo 



^''Hu 



is the standard Fresnel integral. The primes over (^ indicate differentiation evaluated at the 
carrier frequency o)oI27t. 

To study the leading edge of the signal, we again choose To to be effectively infinite. Thus, 
(60) becomes, for the leading edge of the field. 



whereas, for the trailing edge 



2 27r'I^J-^ 



e^[ip"'l2yi^y 



{TTipy/^ 



v{y)dy 



(61) 



L2 27r^/^J-x 



T,-e-(^"'l2Yl'^'y 



{TTifT' 



"y/2 



v{y)dy 



(62) 



As indicated by (61) and (62), the leading edge and the trailing edge of the signal envelope have a 
different form. This asymmetry is a result of the <^'" term in the argument of the Fresnel integrals. 
To characterize completely the envelope function, we need to introduce another parameter /3 
defined by 



;8 = [^"'(a)o)/2]>/-'[7r"(a)o)]- 



1/2 



(63) 



in addition to the "buildup" time th. Gershman [1952] has introduced the fS parameter in con- 
nection with a rather qualitative discussion of reflection of an rf pulse from the ionosphere near 
its critical frequency ojc- For example, if coo/cor/^ 0.9, he estimates that (3 is from 0.23 to 0.3 for 
the E layer and from 0.14 to 0.21 for the F layer. Thus, it is of considerable importance to know 
if the signal envelope shown in figure 2, for /3 = 0, is significantly modified for finite ^ values of 
this order. In an example, Gershman shows that near the steepest portion of the leading and 
trailing edges of the signal the envelope is not modified by more than 2 percent over that shown 
in figure 2. However, differences are somewhat greater in the "forerunner" or "tail" regions. 
A proper understanding of this phenomenon awaits a detailed numerical study of the integral in (60). 
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An interesting limiting case which may be studied without difficulty is when (mp")^/'^ is effec- 
tively zero (which occurs at an extreme of the group velocity). In this case, it is not difficuh to 
show that (61), for (p"'{(0{)) < 0, reduces to 



£"(/:) ^e'^^O^-McoQ) f 



i/.r 



(64) 



while, for ip"'{(0{)) > 0, 



£'(^)=e'^o^-Ma.o)n-r 



2 \i/3 



.^ 



(65) 



where 



n 



1 r-f 



(yjdy 



(66) 



As usual, d is the time measured from the instant —{Tol2)-h(p'{a)o) which, in the absence of dis- 
persion, is the leading edge of the signal. Thus, in (64), the function T characterizes the shape 
of the envelope of the signal when (p'" is negative (i.e., group velocity has a minimum at coo). The 
function T{x) is shown plotted in figure 1 for a range of x from — 5 to -h 3 where the abscissa is 
proportional to time from left to right. On the other hand, in (65), the function I — T characterizes 
the shape of the signal envelope when if'" is positive (i.e., group velocity has a maximum at coo). 
However, in this case, the abscissa in figure 1 corresponds to time moving from right to left. 
Asymptotic expansions for the function T{x) may be readily obtained by making use of formulas 
quoted by Antosiewicz [1964| who also gives some additional numerical data of closely related 
functions. Retaining only the leading term in the asymptotic representation, it follows that 



T{x) ~ 1 



1 



2^1/2^3/4 



exp 



(-!'■'■) 



for x>>l 



T(x) 



^l/2(_^)3/4 



COS {-::(— x)^l^-^—] for- 



r > > 1. 



This points up the marked asymmetry of the envelope function about the origin. As also indicated 
in figure 1, T{x) is highly oscillatory for negative x (i.e., in the "forerunner" region) whereas it is 
monotonic for positive x (i.e., in the "posterior" region). 

The shape of the signal envelope at the traihng edge lor the case ip"{(x)o) = is found in a very 
similar manner to the above. Thus, for (p'"(coo) < 0, (62) becomes 



£'(^)=^/a.o^-Mo>o) J 1—7^ 



while, for (^'"(coo) > 0, 



/r(^) = ^/o;,)^-'V(a>o) 7 



0-To 



e-To 



2 



-^ 



1/3 



^ 



1/3 



(67) 



(68) 



The assumption in the above of effectively infinitely long pulse length requires here that 
To I 2l(p"' I ^''^ > > 1. In this case, it is necessary to suitably superimpose the responses in (64), 
(65), (67), and (68). An interesting consequence of the superposition is that the response or signal 
envelope is not symmetrical about the center of the pulse. This is to be contrasted with the 
situation exemplified by (57) which leads to a symmetrical envelope. In fact, in general, it may 
be asserted that distortion resulting from even-ordered derivatives of (^(w) leads to symmetrical 
pulse distortion. Asymmetry is introduced when any of the odd-ordered derivatives of ip{o)) 
contribute to the final response. 
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8. The Forerunner 

In the preceding discussion of quasi-monochromatic pulse propagation, it is tacitly assumed 
that the spectrum of the signal has a predominant peak at the carrier frequency coo/27r. This is 
usually justified for the main body of the signal provided the dispersion is well behaved in the 
vicinity of coo (i.e., no strong absorption bands near (Oo). However, in any dispersive channel such 
as a waveguide, some energy always arrives immediately after the instant t=dlc where d is the 
distance between the source and the observer and c is the velocity of light. While this part of 
the signal usually has a very small energy content compared with the main signal occurring in the 
region at t = (p'i(Oo), it is of interest to make a quantitative estimate of this "precursor" as it is 
sometimes called. This is readily accomplished by applying the stationary phase method directly 
to (42) where we regard gio)) to be slowly varying compared with other factors in the integrand. 
Obviously, this will be vaUd only at very small times in the response where the very-high frequencies 
dominate. The integral to evaluate these has the same form as (4) and the result (10) may be 
used directly. For example, if the source Eo{t) is given by 

Eo(^) — sin coot for ^ > 

= for ^ < (69) 



then, 



g{(o) = I sin o)ote-'^^dt = -f ^^^^^ , (70) 



COq— co^' 



and the "precursor" is given by 



Eit) ^ (-T i 2^' 21 I ^y\|i/2 ^^^ Mco.)-co.^±7r/4], (71) 



where + 7r/4 is to be used for ip" > and — 7r/4 for ip" < 0, and where co.s is the saddle point deter- 
mined by 

(p'{o)s) = 0. 

Expressing this result in terms of waveguide mode parameters. 



E(t) = \-] —I ^ ^1/ . X COSOJs 

\7r/ co^-cof , , al/i;g(co) . 



\ Vp{o}s)l 4. 



(72) 



where — 7r/4 is to be used for dil/vg)/da) > 0, and -h7r/4 for d(llvg)ldo) < 0. It is understood that 
this result for E{t) is to be used only to estimate the early forerunner response. That is, in the time 
region following t = djc seconds but before the arrival of the main body of the signal near t = ip'{cx)o) 
= d/vgicoo). Because o)s > coo in this early forerunner region, (72) is usually of negligible magnitude 
compared with the field amplitude in the main body of the signal. 

9. Concluding Remarks 

In this paper, an attempt has been made to collect together various methods for handling 
transient calculations for a dispersive channel. The approach has been to employ simphfications 
at the expense of rigor. Like most asymptotic results, the indicated formulas should be applied 
to a given problem only when accompanied with a certain amount of caution. 

1400 



I thank Harold Dougherty for carefully checking all the equations, and Mrs. Eileen Brackett 
for help in preparing the manuscript. The work was supported, in part, by the Air Force Cam- 
bridge Research Laboratories, PRO-62-201. 

10. References 

Antosiewic'z, H. A. (June 1964), Handbook of Vlathcmatical Kunctions, ed. M. Abraniowitz and I. Stegun, NBS Appl. Math. 

Series AMS 55, 435-479. 
Brillouin. L. (1914), About the propagation of Ught in dispersive media, Ann. Physik ser. 4 44,203. 
Budden, K. G. (Jan. 1951), The propagation of a radio atmospheric, Phil. Mag. 42, 1. 

Campbell, G., and R. M. Foster (1948), Fourier Integrals for Practical Application (Van Nostrand Co. Ltd., New York, N.Y.). 
Elhott, R. S. (Oct. 1957), Pulse waveform degradation due to dispersion in waveguide, IRE Trans. Microwave Theory 

Tech. MTT-5, 254-257. 
E«wing, M., W. S. Jardetsky, and F. Press (1957), Elastic Waves in Layered Media (McGraw-Hill Book Co., Inc., New York, 

N.Y.). 
Gajewski, R. (Jan. 1958), Influence of wall losses on pulse propagation in waveguides, J. Appl. Phys. 29, 22-24. 
Gershman, B. N. (1952), On deformation of electromagnetic pulses propagated in the ionosphere, Zh. Tekhn. Fiz. (USSR) 

22, 101-104. 
Ginzburg, V. L. (1964), Propagation of Electromagnetic Waves in Plasma (Pergamon Press, Oxford). 
Ito, M. (Oct. 1964), Dispersion of a step modulated carrier wave in a waveguide, Proc. IEEE 52, 1250. 
Karbowiak, A. E. (Feb. 1957), Propagation of transients in waveguides, lEE Monograph No. 224R. 

Kotvun, A. A. (May 1958), Transient processes in waveguides. Radio Eng. and Electron. (USSR) (English Transl.) 3, 103-122. 
McLachlan, N. W., and P. Humbert (1941), Formulaire pourle calcul symbolique, Memorial des Sciences Mathematiques, 

Fasicule 100, p. 34. 
Pearson, J. D. (Aug. 1953), The transient motion of sound waves in tubes, Quart. J. Mech. Appl. Math. 6, pt. 3, 313-335. 
Proud, J. M.. P. Tamarkin. and E. T. Kornhauser (Jan. 1956), Propagation of sound pulses in a dispersive medium, J. Acoust. 

Soc. Am. 28, 80-85. 
Sommerfeld, A. (1914), About the propagation of light in dispersive media, Ann. Physik ser. 4 44, 177. 
Wait, J. R. (Sept. 1958), Propagation of VLF pulses to great distances, J. Res. NBS 61, 187-203. 
Wait, J. R., and K. P. Spies (Dec. 30, 1964), Characteristics of the earth-ionosphere waveguide for VLF radio waves, NBS 

Tech. Note No. 300. 
Walther, K. (May 1961), Pulse compression in an acoustic waveguide, J. Acoust. Soc. Am. 33, 681-686. 



1 1 . Additional Related References 

Dougherty, H. T., and M. E. Johnson (Sept. 18, 1964), Tabulation of Airy functions, NBS Tech. Note No. 228. 

Gajewski, R. (1956), On transient radiation of a dipole inside a waveguide. Acta Phys. Polon. 15, 25-41. 

Knop, C. M. (July 1964), Pulsed electromagnetic wave propagation in dispersive media, IEEE Trans. Ant. Prop. AP-12, 
494-496. 

Knop, C. M., and G. I. Cohn (Sept. 1963), Comments on "Pulse waveform degradation due to dispersion in waveguide," 
IEEE Trans. Microwave Theory Tech. MTT-11, No. 5, 445-447. 

Rubinowicz, A. (1950), Propagation of a cut-off train of De Broglie waves. Acta Phys. Polon. 10, 79-96. 

Rubinowicz, A. (1954), Uber die Fortplanzung unstetiger elektromagnetischer Signale in Wellen leitern. Acta Phys. Polon. 
13, 115-133. 

Saxton, W. A., and H. J. Schmitt (1963), Transients in a large waveguide, Proc. IEEE 51, 405. 

Schmitt, H. J. (1964), Plasma diagnostics with short electromagnetic pulses, IEEE Trans. Nucl. Sci. NS-11, 125. 

Schmitt, H. J. (Jan. 1965), Pulse dispersion in a gyrotropic plasma, unpublished report from Sperry Rand Res. Center, Sud- 
bury, Mass. 



(Paper 69D1 1-572) 



1401 



